notebook.community

Edit and run



In [1]:

    
%matplotlib inline
from matplotlib import pyplot as plt
import numpy as np



In [2]:

    
fname = 'proj_bias_hist_idx_000_box_000.npy'



In [3]:

    
im_list = np.load(fname)
print im_list.shape









    



(2025, 16, 16, 1)

for idx, im in enumerate(im_list): #if idx > 0: # break implt = plt.imshow(np.log10(np.abs(im[:, :, 0])+1e-6), vmin = -6, vmax = 6) plt.colorbar(implt) plt.show()



In [4]:

    
fname2 = 'full_proj_box_000.npy'#'full_bias_hist_idx_000_box_000.npy'



In [5]:

    
particle_box = np.load(fname2)



In [6]:

    
particle_box.shape









    Out[6]:





(720, 720, 1)



In [7]:

    
np.max(particle_box), np.min(particle_box)









    Out[7]:





(9936.0, 0.0)



In [11]:

    
particle_box.shape









    Out[11]:





(720, 720, 1)

for idx, im in enumerate(particle_box): #if idx > 0: # break implt = plt.imshow(np.log10(im[:, :, 0]+1), vmin = 0, vmax = 4) plt.colorbar(implt) plt.show()



In [9]:

    
from matplotlib.colors import LogNorm
#norm = LogNorm(vmin = particle_box.min(), vmax = particle_box.max())



In [10]:

    
img = plt.imshow(particle_box[:,:,0] - particle_box.min() + 1e-6, norm = LogNorm(vmin = 1e-6, vmax = particle_box.max()-particle_box.min()))
plt.colorbar(img)









    Out[10]:





<matplotlib.colorbar.Colorbar at 0x7f09dc3a3290>



In [ ]: